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Abstract 



We present a coherent theoretical study of uhrarelativistic heavy-ion data obtained at the 
CERN/SPS by the NA35/NA49 Collaborations using 3+1-dinicnsional relativistic hydrodynam- 
ics. We find excellent agreement with the rapidity spectra of negative hadrons and protons and 
with the correlation measurements in two experiments: 5 -I- 5 at 200 AGeV and Pb + Pb at 
160 AGeV (preliminary results). Within our model this implies that for Pb + Pb {S + S) a 
quark-gluon-plasma of initial volume 174 fnr' (24 fnt^) with a lifetime 3.4 fm/c (1.5 fm/c) 
was formed. It is found that the Bose-Einstein correlation measurements do not determine the 
maximal effective radii of the hadron sources because of the large contributions from resonance 
decay at small momenta. Also within this study we present an NA49 acceptance corrected 
two-pion Bose-Einstein correlation function in the invariant variable, Qinv 
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1 Introduction 



In the last decade many experiments have been performed in the attempt to find evidence for the 
existence of a quark-gluon plasma (QGP). Ever higher energies and/or masses have been involved 
in order to increase the lifetime of the system either by increasing the initial energy density or 
the size of the system. The probability of preparing a strongly interacting system that shows 
thermodynamical behaviour and therefore is treatable by well known thermodynamical or fluid 
dynamical methods increases with the size of the system. In this paper we describe a comprehensive 
hydrodynamical study of data taken at the CERN/SPS by the NA35 and the NA49 Collaborations. 
Our investigation takes into account all available data spectra and correlation functions for different 
particle species. We present results of the analysis [|l[] of Pb + Pb at 160 AGeV using relativistic 
fluid dynamics and assuming an equation-of-state containing a phase transition. These results are 
discussed by comparison to our previous flndings for S" + at 200 AGeV which were obtained by 
applying the same computer code HYL ANDER . 

Many hydrodynamical models js], ^, ^ are available which describe the dynamics of relativistic 
heavy-ion collisions. HYLANDER belongs to the class of models which apply S+l-dimensional 
relativistic one-fluid-dynamics. It provides fully three-dimensional solutions of the hydrodynamical 
relativistic Euler-equations p. HYLANDER has been successfully applied at SPS energies to the 
reaction Au + O and especially to -|- S, a reaction to which we refer several times in this paper. 
Here HYLANDER was used to reproduce Q simultaneously mesonic and baryonic rapidity and 
transverse momentum spectra of the S+S reaction at 200 AGeV. Corresponding measurements had 
been performed by the NA35 Collaboration Q. Based on the successful description of the measured 
single-inclusive spectra, predictions were made for Bose-Einstein correlation (BEC) functions |^ IC]. 



Those predictions agree quantitatively with the measurements |11, Q {of. also Fig. 5a). The 
model also reproduces the photon data for S + Au collisions at SPS energies |13| ] and gives a simple 
explanation for the "soft-p_L puzzle" |14] and the complex behaviour of the radii extracted from 



pion and kaon correlations and explains the difference in the extracted radii for pions and kaons in 
terms of a cloud of pions due to the decay of resonances which surrounds the fireball (pion halo) 

In the present paper we exhibit the space-time geometries of the real hadron sources, i.e., the 
freeze-out hypersurfaces. Results ioi S + S at 200 AGeV and Pb + Pb at 160 AGeV are directly 



compared. For the first time we show our results of effective radii compared to the BEC data of the 
NA35 Collaboration. Furthermore, we present a calculation of the detector-acceptance-corrected 
two-pion correlation function C2{Qinv) in the invariant variable, Qmv In doing so we would like 



to bring the readers attention also to refs. 17] where all the features of these specific types of 
two-particle BEC functions were extensively discussed. The following discussion is a self-consistent 
description of two different heavy-ion experiments. Here we reproduce simultaneously single in- 
clusive spectra of negative hadrons and protons and Bose-Einstein correlations of identical pions 
using only one and the same equation of state. The contents of this paper represents a summary 
of many contributions 0, |, |, g [l^, H in an effort to give a description of heavy ion data taken 
at the CERN/SPS in terms of 3+1 dimensional relativistic hydrodynamics. 



2 Modeling the hydrodynamical solutions 



HYLANDER can use initial conditions ranging between the extremes defined by the Landau |19] 



and the Bjorken |2C] initial conditions. One has to specify an equation of state and a set of param- 
eters which describe the initial conditions. In view of our former results for 5 + 5" at 200 AGeV, 
the same equation of state (EOS) is used for the treatment of Pb + Pb at 160 AGeV. It exhibits 
a phase transition to a quark-gluon plasma at a critical temperature Tq = 200 MeV {cf. refs. 

[2^ , |23|]). This EOS, which has no dependence on the baryon density, is plotted in Fig. 1. In 
the following we sketch the basic features of the model which was introduced in ref. Q. Our model 
uses the five parameters, Kl, A, y/\, and cr, which are explained below. 



A reaction of two baryonic fluids leads to a deceleration of the projectile and target baryonic 
currents and thus to the spread of their width in momentum space. For the initial baryon density 
distribution in rapidity we write 



dh _ 



(1) 



where iym give the positions in rapidity of the two maxima and a gives the width of the baryonic 
density distribution after the collision, respectively. Cy is a normalization constant, the value of 
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which can be determined from the requirement that when eq. @) is integrated over the whole 
accessible rapidity interval, it is equal to the total baryon number of the system. A spatial baryon 
distribution (c/. Fig. 2a) is then easily derived by evaluating 

BJz) = _}_^^ (2) 

where R is the mean radius of the initially radially smeared out cylindrical fireball (c/. ref. [jlO|). 
In our model we impose an initial rapidity field y{z) on the fluid. Its modulus is presumed to be a 
function only of the longitudinal coordinate z with its shape constrained by two boundary values: 
The rapidity should vanish at z = and asymptotically reach its maximum value ycm at z = zizta 
[ycm being the center of mass rapidity and ta the time it takes for two nuclei at the speed of light 
to penetrate each other, respectively). We parameterize the function y{z) by the slope parameter 
dy as 

y{z) = ycm tanh [ay \z\] . (3) 

Rather than using ay as a free parameter in our model, we use the two parameters A and t/a 
determining the slope parameter ay in eqn. (^) through 

y{z = ±A/2) = yA • (4) 

In eq. (Q) the quantity A is the spatial longitudinal extension of the initial fireball (in ref. Q called 
the "Landau volume") and y^ is the absolute value of rapidity at z = it A/2. 



The kinetic energy of the two incoming baryonic fluids is converted into internal excitation (ther- 
mal energy) of a third fluid which is created in the central region. The relative fraction Kl of the 
thermal energy inside the initial fireball volume is another free parameter which fixes the initial 
state of the formed fireball. 



In table 1 we show the choice of the five parameters which leads to the reproduction of the experi- 
mental 5-1-5' data taken by the NA35 Collaboration. Additionally, it was there assumed that due 
to experimental uncertainty for the centrality of the collision, only 85% of the total available energy 
and the total baryon number have been observed. In Fig. 2a the initial longitudinal distributions 



of energy density, as well as the rapidity, normalized to their maximum values are plotted against 
the longitudinal coordinate z (since we deal with a symmetric system, only the distributions for 
z > are shown) . Due to the constraints of energy conservation the choice of the initial parameters 
A, yA and a leads to a limited two-parameter space for the variables Ki^ and y^i respectively. The 
limited two-parameter space is shown for S" + in Fig. 2b. The crossed lines indicate the choice of 
the parameter pair for the reproduction of the heavy-ion data. 

For Ph + Ph at 160 AGeV we reduce our five-parameter space to a two-parameter space based on 
the results we obtained for S' + S' at 200 AGeV. We assume that the parameter for the longitudinal 
extension of the central fireball scales with 2R/j compared to the parameter for 5" + 5 at 200 
AGeV 0]; 7 is the Lorentz contraction factor. By coincidence there is a factor two increase in the 
longitudinal extension. A, going from the reaction S + S at 200 AGeV to Pb + Pb at 160 AGeV. 
Also the time, ta, it takes for two nuclei at the speed of light to penetrate each other is increased 
by a factor two. For S + S the initial longitudinal velocity field v^^{z) = ta,nh{y{z)) increases as a 
function of z within the initial fireball almost linearly as in the Bjorken initial condition scenario 
(there VBjoTken{z) = zIt\ r is a time scale which can be identified with to). Since A and ta each 
scale with a factor two, by using the Bjorken scaling argument for the initial longitudinal velocity 
field, the increase of the initial longitudinal rapidity field with the coordinate z is reduced almost 
by a factor 0.5 for Pb + Pb compared to S" + S which results in the same choice for the initial 
parameter = y(-z = for both reactions. 

We have checked that the variation of the width a of the initial baryon density distribution to a 
large extent does not affect the calculation of rapidity spectra. Therefore, the value of o remained 
unchanged. The values of A, y/\ and a are given in table 1. Thus we are left with only two 
parameters: the relative fraction of the thermal energy Kl in the central fireball and the rapidity 
ym at the maximum of the initial baryon distribution. Fig. 2b shows the corresponding limited 
two-parameter space for Pb + Pb at 160 AGeV due to energy conservation. The crossed lines 
indicate our choice for the two parameters which are determined by fitting the rapidity spectrum of 
negative hadrons of the preliminary NA49 data. All other spectra and correlation functions, which 
will be discussed below, are predictions of the model. Due to the preliminary state of the NA49 
data we decided to perform a calculation for the choice of an impact parameter bimp — 0. Fig. 2a 
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shows the initial longitudinal distributions of energy density as well as the rapidity, normalized to 
their maximum values (c/. table 1) as a function of the longitudinal coordinate z. 



3 Discussion of the hydrodynamical solutions 

Once the initial conditions and the equation of state are specified, one obtains an unambiguous 
solution from the hydrodynamical relativistic Euler-equations. In our calculations we assume that 
hadronization occurs for all particle species at the same fixed freeze-out energy density e/. Since 
our equation of state is not dependent on the baryon density, the freeze-out energy density easily 
translates into a fixed freeze-out temperature Tf. Our choice for the freeze-out temperature is 
Tf = 139Mey. 

The choice of a fixed freeze-out temperature, Tj, determines the final space-time geometry of the 
hydrodynamically expanding fireball. In Fig. 3a we have plotted the freeze-out regions at different 
times for directly produced hadrons in the z — r plane for 5 + 5 at 200 AGeV and for Pb + Pb at 
160 AGeV, respectively. Each line represents the freeze-out hypersurface at fixed times t. One can 
see very easily that the solutions we obtained from our numerical analysis represent an evolution 
of initially disk-shaped fireballs which emit hadrons from the very beginning of their formation. 
While the relativistic fluids expand in longitudinal and in transverse directions, the longitudinal 
positions of the freeze-out points increase their distance relative to the center. Due to the effect 
of transverse inwardly moving rarefaction waves the transverse freeze-out positions move towards 
the center of the flreball (c/. also ref. In the late stage of the hydrodynamical expansion the 

hadron-emitting fireballs separate into two parts while cooling down until they cease to emit. 

In Fig. 3b we give the full three-dimensional views of the freeze-out hypersurfaces of directly emit- 
ted hadrons for S'-l-5' at 200 AGeV as well as for Pb + Pb at 160 AGeV. The freeze-out space-time 
geometries show a similar behaviour for both heavy-ion reactions although the longitudinal and 
transverse sizes of the systems differ approximately by a factor two. 
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Above we have argued that the initial parameters A and t/a have geometrical scaling features, while 
a remains unchanged. The only parameters which were chosen freely (except for the constraint 
of energy conservation) are Kl and ym {cf. Fig. 2b). Consistent with expectation the degree of 
stopping and thermalization is higher in Pb + Pb and the amount of thermal energy (represented 
by the parameter Kl) in the central fireball increases from 43% {S + S) to 65% {Pb + Pb). The 
location of the maximum density of the two baryon currents ym in rapidity space is significantly 
shifted into the central rapidity region. The baryons for Pb + Pb are almost stopped. The resulting 
high baryonic density in the Pb + Pb case of 2.14 fm~^ in the center is three times higher than in 
S + S. In the Pb + Pb case 73% of the baryons are initially located in the central region compared 
to only 49% in the 5 + 5 case {cf. Fig.2a). 

In Fig. 3c we give for both reactions three-dimensional views of corresponding transverse velocities 
at freeze-out. The two systems each show a maximum value for the transverse velocities: for S + S 
at 200 AGeV we obtain a maximum transverse velocity n™"^(S') = 0.43 whereas for Pb + Pb at 
160 AGeV we get Uj^^^{Pb) = 0.61. But the slope of the transverse increase of the transverse 
velocity u± is smaller for the system Pb + Pb compared to the system S + S. The reason for this 
different behaviour has its origin in the different choice of the initial parameters: due to a higher 
initial thermal energy density the internal pressure of the system Pb + Pb is increased compared to 
S + S, resulting in a relatively faster longitudinal and transverse expansion of the relativistic fluid. 
The rarefaction waves in Pb + Pb also move inwardly faster and inhibit the creation of an adequate 
large transverse velocity field at freeze-out, yielding a smaller transverse slope for u± compared to 
S + S. 

From our calculations we find that for Pb + Pb a quark-gluon-plasma with an initial volume of 174 
/m^ was formed while for S + S there is a QGP with an initial volume of 24 fm^. The maximum 
lifetime of the fireball is increased from 6.9 fm/c {S + S) up to 14.5 frn/c for Pb + Pb. We observe 
approximately the same behaviour for the lifetimes of the QGP. Whereas in S + S the lifetime of 
the QGP was short (1.5 fm/c), in the Pb + Pb case a QGP persists for 3.4 fm/c. The latent heat 
of the phase transition is of the order of 1 GeV/fm^ {cf. Fig. 1). Taking into account that the 
initial volume of the Pb + Pb system is approximately eight times larger than in the S + S case, 
we note that the lifetime increases slower than does the volume. This effect is due to the cooling 
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by transverse rarefaction waves which are only sensitive to the difference in the transverse radius 
which differs approximately by a factor two from S" + S" to Ph + Ph. 



4 Single inclusive spectra and Bose-Einstein correlations 

In the following we discuss some of the results for the single- and double-inclusive spectra of mesons 
and baryons. All calculations are based on thermal as well as on chemical equilibrium. In both 
types of spectra we include the effect of resonance decays. The influence of partial coherence jl^ 
will not be considered here. 



The calculation of single particle inclusive spectra with HYLANDER was extensively discussed 
in ref. ||^]. In particular, the momentum distributions are calculated in terms of the generalized 
Cooper-Frye formula (see ref. ||2^), where explicitly a baryon and a strangeness chemical potential 
have to be taken into account. These potentials have to be introduced, because the assumption 
of chemical equilibrium requires zero strangeness and (in general) nonzero baryon density at each 
freeze-out point. In detail one has to solve the following system of equations for each surface point 
of given baryon density b in its rest frame: 

^bini{fiB,fJ-s) = b, 

i 

^Sini{fj.B,l^s) = s = 0. (5) 

i 

The index i enumerates all resonances and their anti-particles. In eq. (|5|) hi, Si and ni{^B,IJ's) 
denote the corresponding baryon number, strangeness and number density of the ith resonance, 
respectively. For S + S 8X 200 AGeV we find average values for the baryonic chemical potential 
{i^b) = 284 MeV and for the strangeness chemical potential {^s) = 44 MeV using Tf = 139 MeV 
for the freeze-out temperature. For Pb + Pb at 160 AGeV we find average values (/U^) = 363 MeF 
and {^s) = 69 MeV using the same freeze-out temperature. As a result, we obtain the chemical 
compositions of hadrons at freeze-out listed in table 2 (c/. also ref. 0). We stress that due to 
our choice of the freeze-out temperature of Tf = 139 MeV , there are large resonance contributions 
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(40% — 50%) to the pionic spectra. 



Figs. 4a and 4b show our model calculations for negative hadrons and protons p for S + S 
at 200 AGeV and Pb + Pb at 160 AGeV, respectively. The spectra of negative hadrons h~ are 
made up by contributions from negatively charged pions and kaons (directly emitted ones plus 
those originating from the decay of resonances). The K~ spectra consist mainly of directly emitted 
kaons; less than 10% come from the K* resonance. In the case of the NA35 data we obtained our fit 
parameters (c/. table 1) from a simultaneous fit to the rapidity and transverse momentum spectra 
for the negative hadrons h~ and the protons p (for the S + S reaction the proton spectra do not 
include contributions from A decays Q while they do in the treatment of the reaction Pb + Pb) . 
By fitting the h~ rapidity spectrum of the NA49 data we also made sure that the calculated proton 
rapidity spectrum reproduces simultaneously values observed in the NA44 experiment [^]. In the 
case of the Pb + Pb data the transverse momentum spectra were not involved in finding a parameter 
set for the fit of experimental data and are therefore predictions. We mention that unlike in the 
presentation of ref. |j^], we now account in our calculations for the asymmetric n/p ratio, which 
is 125/82 for Pb nuclei and which results in an approximately 20% smaller proton spectrum than 
already shown in [Q. The feature of higher stopping in the Pb + Pb collision scenario is nicely 
reflected in the rapidity spectra of the protons. While the rapidity spectrum for protons in the 
S+S case has a minimum at y = 0, for Pb+Pb, the rapidity spectrum of protons shows a maximum. 



The calculation of Bose-Einstein correlations (BEC) was performed using the formalism outlined 
in refs. P, 10, 18| including the decay of resonances. The hadron source is assumed to be fully 



chaotic. Here we present results for pion BEC only. Kaon correlations have been discussed in refs. 
[H, |9|, |l^. In order to extract effective hadron source radii, we fit our results to the Gaussian formQ 
which has been widely used by experimentalists for the presentation of their BEC data (for the 
choice of the variables, cf., e.g. ref. pOt| ): 

1 



C2(A;i, k2 



1 + A exp 



"2(^11-^11 + Iside^side + Qout^outi 



(6) 



It should be emphasized that in the present model A does not represent the effect of coherence, 
^If one performs a fit to a BEC function in more than one dimension, eq. ^ does not represent the most general 
expression, because of the existence of an "out-longitudinal" cross term [E8| . 
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but the momentum-dependent effective reduction of the intercept due to the contributions from 
the decay of long-hved resonances (c/. [||,[^]). 



Figs. 5a and 5b show our calculations for the effective radii Rgide ^-i^d Rout as functions of 
rapidity yx and transverse momentum Kj_ of the pion pair compared to the corresponding NA35 



and preliminary NA49 data |11, 12 1, respectively. In order to make possible a comparison of the 
calculated effective radii with the experimentally obtained ones, we have to account for detector 
acceptances. In the case of + at 200 AGeV the effective radii as a function of i/k have been 
calculated at K± = 200Mey, while the effective radii as a function of K± have been calculated at 
UK = 4.0 — Ucm ~ 1.0 . In the case of Pb + Pb at 160 AGeV the effective radii as a function of 
K± were calculated at yx = 4.5 — ycm ~ 1-6. The effective longitudinal radii as a function of 
K± are evaluated in the longitudinal comoving system (LCMS). All our calculations, which in the 
case of S" + S* have been true predictions, agree surprisingly well with the data. In refs. [0, ^, IC] 
we have shown also the effective radii for both types of heavy-ion collisions for yx = 0.0 and 
K± = 0.0. There the effective radii take even larger values compared to the ones we show here 
when comparing them to the data, because the contributions from resonance decays to BEC take 
their maximum values at low momenta. Therefore, we have strong evidence that one cannot extract 
the maximal effective radii of the hadron sources from the Bose-Einstein correlation data obtained 
at the CERN/SPS since the acceptances of the experiments do not obtain data at yx = 0.0 and 
K± = 0.0. It is important to note that in the case of the pion interferometry, the presence of a 
resonance halo increases the size of the fireball in the central region by factors ~ 2.1 (~ 2.0) in 
longitudinal and ~ 1.4 (~ 1.3) in transverse direction for Pb + Pb {S + S). Values for the maximal 
possible resonance halo-size are given in table 1. 



In Fig. 6 we show a two-particle Bose-Einstein correlation function C2{Qinv) as a function of the 



invariant variable, Qmv This BEC function is defined as follows {cf., also refs. 17, |2^). Let 
P2(^i)^2) be the inclusive two-particle spectrum for two identical pions with momenta ki and k2 
and pi{k) the inclusive single-particle spectrum for a pion with momentum k, respectively: 
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The two-particle BEC for a single pair of two pions with momenta ki and k2 takes the form 

C,(h,k2)) = 'f^'^^'l = 1 + 9'^ ''I . (8) 
Pi(fci) • Pi(^2) Pi{ki) ■ pi{k2) 

The separate contributions c{ki,k2), Pi{ki) and pi{k2) have to be calculated for each pair of pion 
momenta including also the interference contributions from the pions which originate from the 
decay of resonances (c/. refs. |l^ |T7|] ). 



The two-particle Bose-Einstein correlation function C2{Qinv) as a function of the invariant variable, 
is given by 



With = k'^- k!^ we have 



C2{Qi 



1 + 



h{Qi 



hiiQi 



(9) 



hliQinv) = J '^^^ J duJ2S[Qinv " \/^q^] Plih) Plih) 
hiQinv) = duJi duJ2 S [Qinv " yf-q^^Qfi ] cQilM) , 



(10) 



and k^^ represents the 4-momentum of a pion. The integrations have to take into account the par- 
ticular detector acceptance for the experiment under consideration. Here we considered the phase 
space 1.1 < y < 2.1 and bOAIeV < k± < QOOMeV, which covers the detector acceptance of the 
NA49 experiment. 



In the previous section we mentioned that in the case of pion interferometry we have to deal with a 
large fraction ~ 40% to 50% - of the pions originating from resonances. This statement refers to an 
average over all particle momenta. In case of the NA49 detector acceptance the main contributions 
to the two-pion BEC come from thermal vr^ (72.8%) and vr^ contributions from p- (15.6%), uj- 
(8.7%) and r/-decays (2.9%). The effects of resonance decays on the two-particle Bose-Einstein 
correlation function C2{Qinv) are also shown in Fig. 6, where we successively have added contri- 
butions from the specific resonance decay channels to the thermal (direct) pion contributions. The 
more resonances we take into account, the narrower the correlation function becomes because the 
addition of a resonance halo increases the effective source size (c/. refs. |l^). As discussed 
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in our earlier papers, the resonance contributions from the ?7-resonance yield an apparent intercept 
reduction. The intercept, Iq = C2{Qinv = 0), of the BEC takes the value /q ~ 1 + 0.97^ 1.94 (c/. 
also ref . ||l^ ) . Here we do not compare our numerical result to the preliminary data of the NA49 
Collaboration |jl2|, because the current experimental correlation function represents a data sample 
of only 529 events. Furthermore, we believe that there may be problems in the preliminary data 
with the experimental two-track resolution and/or Coulomb overcorrections that are apparent at 
very small values of Qinv Nevertheless the present particular calculation demonstrates the com- 
plexity of BEC functions as they emerge from measurements. 

In general, we expect the observed two-particle BEC to be strongly dependent on detector ac- 
ceptances as well on the particular contributions from the decay of resonances. Bose-Einstein 
correlations are a very complicated observable defined through quantum statistics. These functions 
depend per definition on the choice of momenta under consideration. In general, different detec- 
tor acceptances in different experimental setups should lead for the same heavy-ion collision at a 
fixed reaction energy to different results. Thus, the interpretation of BEC measurements is also 
complicated. The interpretation of extracted inverse widths (effective radii) of experimental BEC 
depends on the specific detector acceptance under consideration. Therefore, rather than interpret 
experimental fitted quantities such as effective source radii, we propose the reader should pay at- 
tention to the model which analyzes the data. From the hydrodynamical treatment we learn that 
the hadron source (the real fireball) is represented through a very complex freeze-out hypersurface 
(c/. Figs. 3a, 3b). The longitudinal and transverse extensions of the fireball change dynamically as 
a function of time, rather than show up in static effective radii. A future study [^0| of the BEC for 
both pions and kaons for two different detector acceptances (which could be done by considering si- 
multaneously the experimental results of the NA49 and the NA44 Collaborations) is in preparation. 



5 Summary 

We have shown that data from two different heavy-ion experiments for single and double inclusive 
cross sections of mesons and baryons can be reproduced with a self-consistent three-dimensional 
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relativistic hydrodynamic description assuming an equation of state with a phase transition to 
a QGP. Our data analysis indicates a stronger stopping and an enhanced transverse flow in the 
case of Pb + Pb collisions compared to S + S collisions at CERN/SPS energies. In particular, the 
preliminary Pb + Pb data can be explained by simple scaling assumptions in the initial conditions 
coming from S+S, although the final distributions do not show these scaling features {e.g., compare 
the final rapidity spectra of protons (c/. Fig. 4a)). 

Bose-Einstein correlation functions for pions have also been calculated. The NA35 and NA49 data 
on interferometry are surprisingly well described. Because the largest contributions to BEC from 
resonance decays are at small particle momenta, the current BEC experiments at CERN do not 
measure the maximal possible interferometry radii for identical pions. We also exhibit the freeze-out 
hyper surf aces; these surfaces rather than the effective radii extracted from experimentally measured 
BEC functions represent the true space-time geometries of the sources. This caveat should always 
be kept in mind when trying to interpret the BEC data. 

The results of this work constitute further evidence that heavy-ion collisions in the SPS region show 
fluid dynamical behaviour and can be described by assuming an equation of state with a phase 
transition from QGP to hadronic matter. This result has to be regarded in line with previous 
hydrodynamical studies of the Marburg group ll, g |, |, in and from other groups Of 



course, the EOS we have chosen might not be the only one which is able to describe the current 
SPS heavy-ion data. Further analysis which allows for different equations of state is in preparation 
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Figure Captions 



Fig. 1 Equation of state with phase-transition to a quark-gluon plasma. Plotted are the speed of 
sound, Cq, and the energy density, e, as a function of temperature, T |23]. The curves are a 
result of a fit [^] to lattice QCD calculations pi]]. 

Fig. 2a Initial distributions of energy and baryon density, as well as the rapidity, normalized to 
their maximum values (c/. table l)and plotted against the longitudinal coordinate z. 

Fig. 2b Boundaries in the two-parameter planes for the relative fraction of thermal energy in 
the central fireball and the rapidity y^a at the maximum of the initial baryon distribution, 
respectively. The crossed lines indicate our particular choices for the calculations. 

Fig. 3a Time contour plots of the freeze-out hypersurfaces in the z — r plane. Each line represents 
the freeze-out hypersurface at a fixed time t. 

Fig. 3b Three-dimensional view of the freeze-out hypersurfaces. 

Fig. 3c Three-dimensional view of the transverse velocity fields at freeze-out. 

Fig. 4a Rapidity spectra for negative hadrons {h^) and protons {p). The data points of the NA35 
data are taken from [^. The data points of the preliminary results of the NA49 Collaboration 
for negative hadrons from central Pb+Pb collisions at 160 AGeV are shown with two different 
markers. Circles (diamonds) stand for the measurements from the VTP2 (MTPC) (c/. refs. 

©Hi). 

Fig. 4b Transverse momentum spectra of negative hadrons (/i^) and protons (p). The data points 
of the NA35 data are taken from For Pb + Pb at 160 AGeV the integrations with respect 
to rapidity y have been performed over the interval \y\ < 3.0. 

Fig. 5a Effective radii extracted from Bose-Einstein correlation functions as a function of the 
rapidity yx of the pair and the transverse average momentum Kj^ of the pair for all pions 



compared to NA35 data [y, |l|]. 



Fig. 5b Effective radii extracted from Bose-Einstein correlation functions as a function of the 
transverse average momentum K± of the pair for all pions compared to preliminary NA49 
data [H]. 
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6 Two-particle Bose-Einstein correlation functions C2{Qinv) of negatively charged pions. The 
contributions from resonance decays are successively added to the correlation function of 
thermal tt" (dotted lines). The resultant correlation function of all 7r~ is given by the solid 
line. The integrations have been performed with respect to the detector acceptance of the 
NA49 experiment (see text). 
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Table Caption 



Table 1 Properties of initial fireball extracted from a hydrodynamical analysis of the S + S NA35 



data [| and the Pb + Pb NA49 data ^ 



Table 2 Chemical compositions of mesons, baryons and anti-baryons at freeze-out for S" + S" at 
200 AGeV and Pb + Pb at 160 AGeV . Listed are their relative fractions. 
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S+S 


Pb+Pb 


Fit parameters 


Rel. fraction, Kl, of thermal energy in the central fireball 


0.43 


0.65 


Longitudinal extension, A, of central fireball 


0.6 fm 


1.2 fm 


Rapidity, yA, at edge of central fireball 


0.9 


0.9 


Rapidity, y^, at maximum of initial baryon distribution 


0.82 


0.60 


Width, (7, of initial baryon y-distribution 


0.4 


0.4 


Output 


Center of mass rapidity, ?/cm 


3.03 


2.92 


Max. initial energy density, emax 


13.0 GeV/fm^ 


20.4 GeV/fm^ 


Max. initial baryon density, B'^^^ 


2.66 


3.20 fm-^ 


Rel. fraction of baryons in central fireball 


0.49 


0.73 


Max. lifetime 


6.9 fm/c 


13.5 fm/c 


Initial volume of QGP 


24 frri^ 


174 fm^ 


Lifetime of QGP 


1.5 fm/c 


3.4 fm/c 


Max. halo-size in "longitudinal" direction 


2.5 fm 


5.6 fm 


Max. halo-size in "side" direction 


1.0 fm 


2.4 fm 


Max. halo-size in "out" direction 


1.3 fm 


2.6 fm 



Table 1 
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Particle species 


Mass [GeV/c^] 


s+s 


Pb+Pb 


TT (stable) 


0.139 


51.094 % 


49.201 % 


u 


0.783 


2.452 % 


2.304 % 


ri 


0.958 


0.297 % 


0.279 % 


V 


0.549 


2.941 % 


2.765 % 


P 


0.770 


7.429 % 


6.980 % 


K+ (stable) 


0.494 


7.689 % 


7.763 % 




0.498 


2.573 % 


2.582 % 


A 


1.116 


1.174 % 


1.356 % 


S 


1.193 


2.232 % 


2.578 % 


A 


1.232 


4.115 % 


5.014 % 


N (stable) 


0.939 


9.009 % 


10.960 % 




0.893 


0.952 Vi 


0.951 % 




1.386 


0.613 % 


0.707 % 




1.320 


0.497 % 


0.541 % 


R- (stable) 


0.494 


4.400 % 


3.881 % 


anti-iv:° 


0.498 


1.472 % 


1.289 % 


anti-A 


1.116 


0.049 % 


0.036 % 


anti-E 


1.193 


0.092 % 


0.068 % 


anti-A 


1.232 


0.098 % 


0.066 % 


anti-iV (stable) 


0.939 


0.215 % 


0.147 % 


anti-i^* 


0.893 


0.546 % 


0.479 % 


anti-E* 


1.386 


0.025 % 


0.019 % 


anti-S 


1.320 


0.036 % 


0.028 % 



Table 2 
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